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Recent study on doping effects in the heavy fermion superconductor CeCoIns has shown that a 
small amount of doping induces unexpectedly large broadening of the transition into the high field 
and low temperature (HFLT) phase of this material. To resolve this observation, effects of quenched 
disorder on the second order transition into a longitudinal Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state are examined. The large broadening of the transition is naturally explained as a consequence 
of softness of each FFLO nodal plane. The present results strongly support the scenario identifying 
the HFLT phase of CeCoIns with a longitudinal FFLO state. 
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Understanding the high field and low temperature 
(HFLT) phase [l|-[3| realized in the heavy-fermion super- 
conductor CeCoIns in the parallel field configuration is 
an issue under hot debate on unconventional supercon- 
ductivity. Based on a large paramagnetic depairing of 
this material and on a close examination on the phase 
diagram, the HFLT phase has been identified [E with 
a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) (3,13 vortex 
lattice with a modulation parallel to the field [6], which 
will be dubbed hereafter the longitudinal FFLO state. 
This FFLO scenario naturally explains most of proper- 
ties relevant to the HFLT phase, such as the striking 
reduction of the vortex tilt modulus in the HFLT phase 
0, 0] and the second order character of the transition 
(HAT) between the HFLT and the ordinary Abrikosov 
lattice Q phases. On the other hand, following a recent 
observation of a transverse [1[ antiferromagnetic (AFM) 
order, much attention have been paid to magnetic prop- 
erties and have led to an alternative scenario against the 
FFLO one, in which HAT is identified with the onset 
of a Q-phase, which is a pair density wave with no lon- 
gitudinal modulation and accompanied by a bulk AFM 
order [j| [l(| (See also Ref.f f). It has been speculated 
[l2j and demonstrated [l3j], however, that the observed 
AFM order may be an event localized around the trans- 
verse nodal planes in the longitudinal FFLO state [l^ . 
Besides this AFM ordering, a striking doping effect on 
HAT has been reported in heat capacity measurements 
d: For both of AFM dopants (Hg and Cd) and a 
nonmagnetic (Sn) one, even an extremely small amount 
of doping has induced a transition broadening and a dra- 
matic reduction of the heat capacity jump at HAT, sug- 
gesting that an ordering occurring through HAT is highly 
fragile. This is qualitatively and quantitatively different 
from the familiar electronic impurity effects on a uncon- 
ventional pairing transition [l6| and a FFLO transition 
[l7| in which just a shift of the transition point is ex- 
pected. 

In this work, we examine quenched disorder effects on 
HAT between the longitudinal FFLO and the ordinary 
Abrikosov vortex lattices and explain the observed broad- 
ening of the heat capacity near and below HAT 0, [H| 



with the help of numerical details of a microscopically 
derived Ginzburg-Landau (GL) Hamiltonian 7]. A soft 
tilt rigidity 0] of the FFLO nodal plane is found to be 
the main origin of the dramatic broadening of heat ca- 
pacity curves fbH . [Til ] . For comparison, the same analysis 
is also performed for a GL model appropriate for the Q- 
phase scenario [9l4ll1| of the HFLT phase. In the latter, 
the primary effect of quenched disorder is always a simple 
shift of the transition point unaccompanied by a notable 
impurity-induced broadening. Based on these results, 
the validity on the picture (flH, @ identifying the HFLT 
phase with the longitudinal FFLO state is stresseed. 

We consider the Hamiltonian T~L = 
N{0){2ir£, ) 3 T? [h + hp], where 
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are dimcnsionless and valid in the Larkin-Ovchinnikov 
(LO) vortex state with the pair-field (superconducting 
order parameter) A(r) = y/2A(x, y)cos(Qoz + s(r)) near 
HAT. Here, N(0) is the density of states on the Fermi 
surface in the normal state, T c is the zero field supercon- 
ducting transition temperature of the undoped system, 
Q = Qo + d z s is the FFLO order parameter express- 
ing the inverse of the local period of FFLO modulation 
parallel to the applied field H || z, s = sz is the displace- 
ment field of the nodal planes lying in the x-y plane in 
equilibrium, u _L z is the compressional displacement of 
the vortex lattice arising from A(x, y), Hij is an elastic 
matrix of vortices to be defined later, and /?, e > 0. In 
cq. (JIJ , the first two terms describe the mean field order- 
ing of the longitudinal FFLO state, while other terms 
expressing the elasticity of the nodal planes and the vor- 
tices and the coupling between them have been exam- 
ined elsewhere Q- Any length was already normalized 
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by 27r£o, where £o is the coherence length in T = 
limit. In calculating the heat capacity, numerical data 
of the dimensionless coefficients a, 7, and e will be 
used (see below) [l8| . Among possible roles of impurities, 
we focus hereafter on their quenched disorder effects on 
the order parameter fields described by eq.(2), and the 
electronic impurity effects will be commented on at the 
end of this paper. The random field terms in eq.(2), for 
instance, follow from the conventional random T c term 
~ / w(r)|A(r)| 2 in the superconducting GL Hamilto- 
nian Q , and the presence of the factor Qo in the nodal 
plane pinning term proportional to h(r) may be justified 
from the model w(r) = ^ro^^fr — r v ). A possible 
randomness of a implying spatial inhomogenuities of the 
HAT temperature is of a higher order compared with the 
/i-term and thus, was neglected. 

To examine the free energy density T in the FFLO 
state near HAT, thermal fluctuations of s and u will be 
neglected. Then, the method used by Larkin and Ovchin- 
nikov [l9] for a second order transition will be adopted 
and extended here in a self consistent manner to obtain 
T . After taking variations of H with respect to Qo, u, 
and s and keeping the contributions up to 0(s 2 ) in the 
Q 4 term, the variational equations 



aQ + PQl+[h+j(Il- 1 ) ij d z difj)s = -Q 



3/3(5 zS ) 
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and 



[(a + 3(5Q 2 )d 2 z + QiYi^hd'Ad, + edi]s 



(4) 



follow. By taking the random average of eq.©, the self 
consistent equation determining Qo 



PQl = 



h 



a 



1 + 3/2 



(5) 



is obtained, where 



h = 
h = 



d 3 k S 

d 3 k 
WT 3 
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and D k = [I 1 +P(2-3I 2 - 1 2 U- 1 (k)kl/(3)Q 2 }k 2 z +ekl. 
Here, Hij(k) was isotropized in the way SijH with n = 
H 2 k 2 /{AnN({))T 2 ), S denotes the random average of 
h\Ji-\s_ + Il~ 2 k 2 \k ■ ffc| 2 /8, and, for brevity, its |k| depen- 
dence has been neglected. After the k-integrals in eq.®, 
an upper momentum cutoff k c of order unity will be ab- 
sorbed into the bare disorder strength via its redefinition, 



8ok c /(4:ir 2 ) —> 5q. The resulting expressions can be sim- 
plified further because the 7-dependent terms appearing 
through the u-variation are quantitatively negligible. In 
fact, using typical values of 7 (~ 0.015) and j3 (~ 0.004) 
obtained elsewhere [g, [D] through the study of the phase 
diagram of CeCoIn 5 near H c2 (Q), 7 2 N (0)4ttT 2 / (6H 2 /3) 
is estimated at most as ~ Vd~ 2 k^,T c j Ep, where Ep is the 
Fermi energy. For this reason, any 7-dependent terms 
will be neglected below. Then, we have 
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and I2 = —dli/dm, where m = I\ + (2 — 3I2)/3Qq. 

It should be stressed that, although effects of the vor- 
tex displacement are negligible, the presence of the field- 
induced vortices is not negligible because the nonvanish- 
ing £ (~ 0.002) used in our analysis is a consequence 
of the orbital pair-breaking. On the other hand, the soft 
nodal plane implied by such a small e is a reflection of 
the large paramagnetic depairing in CeCoIns and, in the 
impure case, has a crucial impact on the thermodynamics 
near HAT. Note also that the nodal plane is softer in more 
anisotropic systems such as the organic superconductors 
because a mass anisotropy enhances the paramagnetic 
effect in the same field configuration. 

Based on the expressions obtained above, the change 
of free energy accompanying the FFLO ordering follows 
simply from the random-average of "H, and J- is given by 



2/3 y ' c (1 + 3/ 2 ) : 



:(a-h) 



(8) 



The resulting heat capacity C(T) curves below T s (5o) are 
given in Fig.l as a function of Sq, where the heat capac- 
ity is normalized by its value at T s in the pure (5q = 0) 
case in the form C(T) = C(T)/C(T = T S ,S = 0), and 
T S (S ) is the HAT temperature. The fact that C(T) is 
significantly broadened and depressed by a small amount 
of disorder seems to be consistent with the feature seen 
around HAT in CeCoIn 5 [3 EE)- According to eq.Q, 
this sensitivity of C to quite a small 5q is a consequence of 
the fact that the effective disorder strength is not Sq but 
S /e. A smaller e enhances the disorder effect and leads 
to a more dramatic broadening of the transition. Fur- 
ther, note that a large depression of C-value also implies 
that the period ~ Qq 1 of the FFLO modulation remains 
macroscopically long even at lower temperatures. 

Two crucial features are seen in this figure in relation 
to the broadening : First, T s increases with So up to 
~ 10 2 SoT c /e, and further, a broad peak appears in C(T) 
which, as far as Sq/s > 2.5 x 10 -5 , lies much below T s . In 
fact, this broad peak occurs in the region where to > e, 
and hence, tranverse spatial variations (|kj_| 3> \k z \) of 
the nodal planes are dominant. Thus, the broadening 
accompanied by the suppression of the peak in Fig.l is a 
consequence of the softness of each nodal plane. We note 
that the present approach is not applicable to the so- 
called critical region in the close vicinity of T s because of 
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FIG. 1: Results of normalized heat capacity C(T) (thick solid 
curves) below T a (8o) and at a fixed magnetic field following 
from eq.© for 5 = 8.5 x ICT 8 (top), 3.0 x f(T 7 , 1.0 x 10 -6 , 
and 2.0 x 10" 6 (bottom). Data in Ref.0] on the coefficients in 
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eq.(l) at H = 0.5^°' (0) have been used, where H££ ; (0) = 
O.560o/(27rfo) is the 2D orbital-limiting field, and </>o is the 
flux quantum. The right end of each curve corresponds to the 
result at each T s (5q). The thin solid line denotes C(T) below 
T a = 0.354T C in the pure (So = 0) case. 



the neglect of nonlinear corrections in eq.(2), which may 
play important roles when m <C So/e, or equivalently 
\T S - T\ < 1.5 x 10~ 2 T C . In fact, other nonperturba- 
tive approaches are necessary to describe the behaviors 
in T > T s including the Griffiths regime 20]. However, 
the disorder-induced broadening, which is our main fo- 
cus, is a feature far below the critical region and thus, 
is believed to be reasonably captured within the present 
approach. 

To explain whether or not the impurity-induced broad- 
ening seen in Fig.l is peculiar to the FFLO to Abrikosov 
lattice transition, it is instructive to repeat the same anal- 
ysis in the Q-phase scenario (9|-[TT| on HAT of CeCoIn5 
in which the HFLT phase is identified with a spin-triplet 
pair-density wave accompanied by an AFM order with a 
transverse modulation wavevector Q. In this case, after 
being minimized with respect to the coupled AFM order, 




FIG. 2: Results corresponding to Fig.l in the Q-phase sce- 
nario. The same a as in Fig.l is used, while 7+ is T- 
independent. Although T s of the S' = 2.0 x 10~ 3 curve is 
almost the same as that of the So = 2.0 x 10~ 6 one in Fig.l, 
no notable broadening is seen in this figure. 



the resulting GL Hamiltonian takes a form typical of a 
"two band" superconductor expressed by the d-wave pair 
field and a triplet pair field \& with no Q-modulation, in 
which the ordered d-wave vortex solid plays roles of a 
periodic potential for 'F. Thus, an energy cost due to 
a mismatch between the vortex positions of the rf-wave 
pair field and *F leads to a gradient term for the latter. 
If we only have to focus on the amplitude fluctuation 
-0 + of <F (see below), the resulting GL Hamiltonian for 
slowly varying components of >F is equivalent to the fa- 
miliar random T c model for the Ising spin system and is 
expressed by 



H u = iY(0)(2< ) 3 T 2 



+ 3p'* 2 )iPl + 7+(W>+) 2 + 2fc , * ^+ 



(9) 



to the leading orders in where ^0 is a real and 
constant amplitude of the equilibrium lattice solution of 
<F, and the length scales were isotropized. The phase 
fluctuation of 'F neglected above is that of "interband" 
Josephson phase which remains massive 22]. It leads 
only to a small T-independent contribution to a in the 
first term of eq. @ , which can safely be omitted in exam- 
ining the heat capacity. Equation © may be obtained 
in the Abrikosov's mean field analysis [8J and its exten- 
sion focusing on the lowest Landau level modes of 
Then, the expression corresponding to eq.© can easily 
be obtained, which becomes 



1 + 3/4 



(10) 



where the integrals I' n (n = 1, 2) are given by I[ — S' (l — 
7r(m'/(7+A; 2 )) 1 / 2 /2) with an upper momentum cutoff k c , 
and I' 2 = —dl'i/dml , where 5q is the random average of 
h' k h'_ k multiplied by A: 3 /(2tt 2 ), and m' = a + 3/3' In 
addition, the free energy density is given by eq.© with 
I n replaced by I' n . The resulting heat capacity curves in 
the fo ^ state are given in Fig. 2. Although, for sim- 
plicity, the combination 7+fc 2 has been set to be unity 
in Fig. 2, qualitative results are independent of the 7+- 
value because it can be absorbed into the effective value 
of disorder strength. The figure shows that no dramatic 
broadening occurs even for a moderate strength of dis- 
order, and a deviation from the pure (5' — 0) case is 
noticable just near the critical region. Clearly, the main 
effect of disorder is to shift the transition point. 

Figures 1 and 2 show that the transition broade ning 
induced by a small amount of dopings in CeCoIns [Til.1151] 
is consistent only with the case in which its HFLT phase 
is a vortex lattice with a longitudinal modulation. As 
already mentioned, the broadening is seen experimentally 
irrespective of the type of the doped element. This fact 
indicates that the broadening should be ascribed not to 
a change of electronic details but rather to a quenched 
disorder effect due to the dopants. 
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On the other hand, the measurements show that de- 
pendences of the nominal position of HAT on the type 
of dopants are nonuniversal : For the magnetic doping 
such as Hg and Cd, the position of the broad peak of 
C(T) is shifted rather to higher temperatures with doping 
(14j , while the doping of the nonmagnetic element Sn has 
shifted the broad peak to lower temperatures [l5j]. Ac- 
cording to Fig.l, however, the peak position lies markedly 
below the actual transition point T s (8q), and this devi- 
ation between those two temperatures is enhanced with 
increasing S . Besides this, the conventional electronic 
impurity effect due to a nonmagnetic doping induces a 
suppression of the transition temperature into an uncon- 
ventional pairing state [l|| or the FFLO state [23[ . In- 
corporating this electronic effect on T s in eq.(l) is easily 
performed simply by redefining a and does not affect the 
transition broadening in Fig.l induced by quenched dis- 
order. Then, due to the coexistence of such two impurity 
effects on T s competitive with each other, it is not easy 
to predict the actual HAT point in real materials. At 
least, the reduction of the broad peak position due to 
the Sn doping [l5[ suggests the presence of a remark- 
able T s -reduction of an electronic origin [23]. For the 
magnetic doping [l4j], however, the broad peak position 
seems to increase with doping, suggestive of the presence 
of an electronic mechanism leading to a slight increase 
of T s . A separate study on electronic details is needed 
to explain the increase of the broad peak position and of 
the corresponding H C 2(T) [l| due to a magnetic doping. 

The present result also has implication on the issue of 



FFLO phases in strongly anisotropic organic supercon- 



ductors [24J. In these materials, a longitudinal FFLO 



state like that in CeCoIns has not been observed, and 
instead, only a transverse modulation in the plane per- 
pendicular to the field has been argued to appear [24j. 
Theoretically, however, the possible FFLO vortex state 
at higher temperatures is predicted to be the longitudinal 
one Q. The present result leads to the conjecture that, 
due to the strong anisotropy, the e- value in those sys- 
tems which are closer to the vortex-free Pauli limit than 
CeCoIns is extremely low so that the transition into the 
longitudinal FFLO state is dramatically broadened and 
becomes invisible due only to a small amount of impu- 
rities. In other words, the period of the spatial modula- 
tion parallel to the field remains significantly long even at 
lower temperatures. Further discussion on organic mate- 
rials based on the present picture will be given elsewhere. 

In conclusion, the doping-induced large broadening of 
the transition between HFLT and the ordinary Abrikosov 
lattice phases in CeCoIns definitely shows that, in con- 
trast to the proposed Q-phase, the HFLT phase has a 
spatial modulation parallel to the field, and that the lon- 
gitudinal FFLO vortex lattice is its best candidate of such 
a highly fragile HFLT phase. The present result is also 
relevant to the issue on FFLO phases to be observed in 
organic materials in high fields. 

This work was supported by Grant-in- Aid for Scientific 
Research [Grants No. 20102008 and No. 21540360] from 
MEXT, Japan. 



[1] A. D. Bianchi et al., Phys. Rev. Lett. 91 187004 (2003). 
[2] T. Watanabe et al., Phys. Rev. B 70, 020506(R) (2004). 
[3] M. Kenzelman et al, Science 321, 1652 (2008). 
[4] P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964). 
[5] A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. JETP 

20, 762 (1965). 
[6] R. Ikeda, Phys. Rev. B 76, 134504 (2007). 
[7] R. Ikeda, Phys. Rev. B 76, 054517 (2007). 
[8] A. A. Abrikosov, Sov. Phys. JETP 5, 1174 (1957). 
[9] D.F. Agterberg, M. Sigrist, and H. Tsunetsugu, Phys. 

Rev. Lett. 102, 207004 (2009). 
[10] A. Aperis et al., J. Phys. Condens.Matter 20, 434235 

(2008) . 

[11] G. Koutroulakis et al, larXiv:0912. 35481 

[12] R. Ikeda, Phys. Rev. Lett. 102, 069703 (2009). 

[13] Y. Yanase and M. Sigrist, J. Phys. Soc. Jpn. 78, 114715 

(2009) . 



[14] Y. Tokiwa et al., Phys. Rev. Lett. 101, 037001 (2008). 
[15] R. Movshovich et al., presented in M 2 S 2009. 
[16] A. A. Abrikosov and L. P. Gor'kov, Sov. Phys. JETP 12, 
1243 (1961). 

[17] A. B. Vorontsov, I. Vekhter, and M. J. Graf, Phys. Rev. 

B 78, 180505 (2008). 
[18] In the notation in Ref. 7], a = c (2) (T c -y 1/2 r B /v F ) 2 , and 

£ = E 2 /2. 

[19] A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. JETP 
34, 651 (1972). 

[20] G. Tarjus and V. Dotsenko, J. Phys. A 35, 1627 (2002). 

[21] G. Eilenberger, Phys. Rev. 164, 628 (1967). 

[22] A. J. Leggett, Prog. Theor. Phys. 36, 901 (1966). 

[23] H. Adachi and R. Ikeda, Phys. Rev. B 68, 184510 (2003). 

[24] S. Uji et al., Phys. Rev. Lett. 97, 157001 (2006). 



